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Abstract. We study the magnetic relaxation of a system of localized spins interacting through weak dipole 
interactions, at a temperature large with respect to the ordering temperature but low with respect to the 
crystal field level splitting. The relaxation results from quantum spin tunneling but is only allowed on sites 
where the dipole field is very small. At low times, the magnetization decrease is proportional to y/i as 
predicted by Prokofiev and Stamp, and at long times the relaxation can be described as an extension of a 
relaxed zone. The results can be directly compared with very recent experimental data on Fes molecular 
clusters. 

PACS. 75.45.+j Macroscopic quantum phenomena in magnetic systems - 75.40.Mg Numerical simulation 
studies - 75.50.Xx Molecular magnets - 61.46.+W Clusters, nanoparticles, and nanocrystalline materials 



1 A novel relaxation mechanism 

The relaxation of a quantity m(t) toward its equilibrium 
value m(oo) is usually described by a standard model 
taken from a surprisingly limited set, in contrast with the 
extreme diversity of the physical problems to which each 
of these models apply. The simplest model is the linear 
equation 



The weak dipole interaction between different molecu- 
lar groups is not sufficient to produce a magnetic order at 
any accessible temperature and, at thermal equilibrium, 
the +s and -s spins are randomly distributed with an av- 
erage value ro(oo) which depends on the external field. In 
a typical relaxation experiment at very low temperature, 
all spins are initially in the state 



= —a[m(t) — m(oo)] 



(1) 



dm{t) 
dt 

where a is a constant. 

The compound [(tacn) 8 Fe808(OH) 8 ] 8+ (where tacn is 
1-4-7-triazacyclononane) pertains, at low temperature T 
< IK, to a new class which has not been studied until very 
recently This material ||, hereafter called 'Fes', is a 
paramagnet and the quantity m(t) of interest is the mag- 
netization per spin, more precisely its component along a 
well defined axis z which is an easy magnetization axis. 
According to current knowledge (^] , this material is made 
of molecular groups, each of which contains 8 ions Fe +++ 
and has at low temperature a spin s = 10 which results 
from a strong exchange coupling between the 8 Fe +++ 
ions. At low temperature each molecular spin Si is ori- 
ented along z, with the value 



so that m(0) = — s. 

Above 1 K, the relaxation is well described by (jl]) 
the relaxation is exponential H , 



m(t) 



moo 



[m(0) — to(oo)] exp(— t/r) 



(3) 



and 



(4) 



However, at low temperature, the relaxation is not ex- 
ponential. At long times t, it is pretty well described by a 
stretched exponential H 



m(t) — m(oo) = [m(0) — m(oo)] exp[— (t/r)^ 1 } (5) 

where the exponent Pi depends on T and takes the value 
1 above 1 K. Below 0.4 K, fli is nearly constant ||] and 
equal to 



Sf = Si - ±s 



(2) 



Pi w 0.4 



(G) 
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In the present work, kinetic Monte-Carlo simulations 
of the magnetic relaxation of Fes are reported; they are 
made on the basis of a model described in the next section. 



2 The incoherent tunneling model 

At the low temperatures of interest, thermal activation is 
impossible and the relaxation takes place by tunneling. 
In weak external field, tunneling takes place between the 
two states (^), and is only possible for a spin Si if the 
z-component Hi of the local magnetic field is very close 
to 0, say 

-Hi< Hi < Hx (7) 

where Hi is a constant which will be precised below. The 
local field is the sum of the external field H ext and an 
internal field. In most of materials, the internal field is 
partly due to nuclear spins (hyperfine interactions). When 
the magnetic particles are Fe, the hyperfine contact in- 
teraction is nearly absent because the most common Fe 
isotope has no spin, the magnetic isotope has a weak con- 
centration (2%) and moreover a weak spin (1/2). In the 
present work, hyperfine interactions will not be explicitly 
taken into account in the evaluation of the internal field. 
They are but implicitly included in the theory since they 
are probably |l],g] responsible for the resonance width H±. 
This width has been measured by Wernsdorfer et al. [Q in 
Feg and is of order 10 Oersteds. 

Thus, the internal field only depends on the molecular 
spins through the formula 



H 



W 

dip 



(8) 



This field will hereafter be called 'dipole field' although 
it may also contain an exchange component. At long dis- 
tance this exchange part vanishes and the coefficients gij = 
g(Ti — Tj) behave according to the formula 



ff( r ) 



K 



(l-3z 2 /r 2 ) 



(9) 



where z is the component of r on the z axis. 

In a finite system, tunneling is an undamped oscilla- 
tion between two states and does not really result in relax- 
ation. For instance, if the spin Si were isolated, it would 
oscillate in zero field between the two states (||). This is 
not true for the large systems which are studied in experi- 
mental physics, and it is reasonable to make the following 
assumption. 

Basic assumption. 

Any spin i subject to a local field Hi has a probability 
rAHi) per unit time of transition between the two states 
(|2j). The real, nonnegative function rj(Hi) is negligible if 
Hi does not satisfy (pj). 

The choice of the function rj{H) is presumably not es- 
sential provided the above properties are satisfied. A pos- 
sible choice is 



rj(H) = 7,(0) 



1 



H 2 



Hi 

r,(H) = (\H\ > H x 



{-H x <H <H X 



(10) 



The above assumption defines the 'incoherent tunnel- 
ing model' which is studied in the next sections. Since 
the spin-flip transitions modify the fields, the relaxation 
defined by the basic assumptions and by formula (Q) is 
difficult. It has been investigated by kinetic Monte Carlo 
simulations by Prokofiev & Stamp jlj and by Ohm & 
Paulsen More detailed results are presented in sec- 
tion |[ They include a detailed analysis of the effect of the 
sample shape and crystal structure. 

Prokofiev and Stamp have given an analytic descrip- 
tion of the short time behavior. A critical summary of this 
theory is given in the next section. 

In higher external field, the basic assumption must 
be modified since tunneling becomes possible from state 
S\ = — s to an eigenstate (or rather, nearly eigenstate) 
S\ — to, with < to < s, if the external field is such 
that two eigenvectors of the spin Hamiltonian have nearly 
the same energy (resonance condition) ||. Then, the spin 
emits phonons and goes to state S\ = s. It would be easy 
to include this possibility in the following calculations, but 
for the sake of simplicity it will be ignored. This is correct 
if the external field is not too large. 



3 The Prokofiev-Stamp approximation for a 
sphere 

Since the dipole interactions (Q) are essential in the model, 
the shape of the sample should be important if the model 
is correct. The simplest case is that of a spherical sample. 
Then (||) implies, at t = 0, that the local field Hi has the 
same value H(0) for all spins i, except near the surface. 

It will be assumed in this section that the sample is 
spherical and that the external field is such that H (0) = 
0. These conditions allow the rapid reversal of a certain 
amount of spins. Because of this process, the local field is 
no longer uniform, so that the reversal of most of the other 
spins is hindered. The distribution P(t; H) of the local 
fields at time t will be assumed to be continuous, and the 
proportion of spins which can reverse in the time unit is 
of order r/(0)HiP(t; 0); this value is exact if r)(H) =const. 
in a range of width Hi and zero outside. If one assumes 
that the local field Hi and the spin Si are independent at 
each time, it follows from Eq. (10) 

lm(t) = -^t ? (0)ff 1 P(t;0)m(t). (11) 

The function P(t; H) is characterized by its width r(t), 
so that (11) can be replaced by 



dt 



m(t) 



-r)(0)H im (t)/r(t) 



(12) 
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In order to obtain a closed evolution equation for m(t), 
one should relate r(t) to the magnetization m(t). When 
to has its saturation value — s, as is the case at t = 0, 
the width r vanishes for a spherical sample. According to 
Prokofiev and Stamp, r(t) = r(m(t)) with 



r(m) = (1 - \m\/s)H d 



(13) 



where the constant Hd has the order of magnitude of the 
maximum dipole field. 

Relation ( |l3|) can be justified as follows. At short times, 
the local field is the sum of the initial field, which has 
been assumed to be 0, and the dipole field produced by 
those spins which are already reversed. Let £ be the aver- 
age distance between those spins. The typical value of the 
resulting dipole field, as given by (|), is of order K/£ 3 w 
Hda 3 /£ 3 , where a is the distance between spins. The width 
r(m) should have the same order of magnitude, r(m) « 
Hdd 3 /£ 3 . Now, a 3 /£ 3 is the proportion of reversed spins, 
i.e. (1 - H/s)/2. Relationfy) follows. 

Insertion of (|l3| ) into jl^ j yields 



[l-\m(t)\/s}-[l-\m(t)\/s]^ V (0)H 1 \m(t)\/(sH d ) (14) 



the solution of which is 



r]{Q)H x t/H d « \m(t)\/s - 1 - \n[\m(t)\/i 



(15) 



For short times, ( |l5| ) reduces to the result of Prokofiev 
and Stamp M 



1- \m{t)\/ S K y/t/Ta 

while, for long times, ( |l5| ) reduces to 

t 



\m(t)\/s ~ exp 



1 



where 



(16) 



(17) 



(18) 



Formula ( [Tq ) will be seen to describe correctly the 
short time behaviour, while the long time behaviour will 
be seen to be described by (||). Formulae ( p"5| ) and (|l7|) will 
be seen not to be satisfactory at long times. This remark 
shows that the approximations made in the derivation of 
( |l6| ) are only valid at short times. 



4 Non-resonant fields and non-spherical 
samples 



at saturation (i.e. if Si = —s) if the external field is absent. 
However, the replacement by an integral is not exact, and 
the local field is not uniform at saturation, although it is 
very sharply peaked (Fig. |l]a). As a matter of fact, it is 
impossible to cut a perfect sphere in a crystalline material, 
and the edges can only be approximately spherical R . 

Relaxation can thus occur near the surface even if the 
resonance condition Hi = is not satisfied in the bulk. 
The relaxation can propagate and possibly become total as 
t — > oo. Whether this occurs or not is one of the questions 
to be answered by the simulations reported below. 

If the sample is not a sphere (nor an ellipsoid) the 
resonance condition can only be satisfied in apart of the 
sample, even at t = 0, when Si — —s (Fig. yb) so that 
the magnetization has a broad distribution (Figs. [l]c,|]a). 
According to Prokofiev and Stamp jjj , the variation of the 
magnetization at short time t is still proportional to \ftfr 
as in (16); such prediction has recently received beautiful 
experimental confirmations |&0. Moreover in Ref. M is 
found that the short time relaxation constant r contains 
the volume where the local field has its resonance value. 
This prediction ignores the possible extension of this vol- 
ume in time, and is therefore somewhat speculative, but 
it opens the possibility to experimentally probe the field 
distribution ||. 

For a non-spherical sample, the square root law is 
observed on a broader field interval since the field dis- 
tribution is broader. Moreover, the square root law can 
be generalized to any value of the initial magnetization 
namely Q 

m(0) - m(t) w Const x y/t . (19) 

This result holds only if the density of the local field 
does not change much for a field variation of the order of 
the resonance width H±. This excludes spherical samples. 

Some of the formulae written in section ^| for a spher- 
ical sample can be generalized. For instance, if one intro- 
duces the local magnetization m(r, t) and the local distri- 
bution P(r, t; H) of internal field, ( Jll| ) can be generalized 
as 

^m(r, t) = -^(OjfliPfo t; 0)m(r, t) . (20) 

Integrating P(r,t;H) on r, one obtains the distribu- 
tion of local fields at time t, 

p(t; H) = ( d 3 r [ dH'P(r, t; H)8{H' — H) (21) 

J J —oo 

which is shown for t — for various sample shapes on 
Figures |, | and |. 



The analytic formulae ( |15| ) to (g_7() have been obtained for 
a spherical sample when the external field is such that, at 
t = 0, the 'resonance condition' Hi — is satisfied for all 
spins i. 

If in (||) the sum is replaced by an integral, the local 
field is uniform in a spherical sample, and vanishes exactly 



5 Monte-Carlo simulations: model and 
method 

The model sketched in section || is not yet completely 
defined. The crystal lattice, for instance, has not been de- 
fined. 
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The material Fes has a complicated, triclinic crystal 
structure. It can be taken into account in the simulations, 
but this does not warrant realism. A complete, realistic 
calculation would imply a calculation of the coefficients 
gij taking into account the exact shape of the electronic 
wave functions. This would be a difficult task. Moreover, 
the presence of short range exchange interactions, which 
would be more difficult yet to calculate, cannot be ex- 
cluded. 

We believe that all these unknown effects can be taken 
into account by a single parameter, which is the value Hq 
of the external field at resonance for a sphere or at the 
center of a cube. For a cubic sample of a cubic crystal, 
Ho = by symmetry. For the real material Fes, Ho does 
not vanish, cannot be calculated for the reasons explained 
above, but can easily be measured. When it is known, 
the field distribution in a sample of any shape can be 
calculated for a cubic lattice, and then obtained for the 
real material by shifting the field scale by an amount Ho- 
The above views are supported by figure || which shows 
the internal field distribution for a cubic lattice compared 
to other possible models. The distributions are not very 
different apart from a shift. 

Most of our simulations have been performed on a cu- 
bic lattice. We have done a single simulation on a system 
with the crystal structure of Fes (but with localized spins) 
and checked that there is no significant difference with 
simulations on a cubic lattice 

An important parameter is the resonance width H\. 
If it were larger than i/dip, the relaxation would be fast 
and exponential. The experimentally observed, slow re- 
laxation is certainly related to a small value of H\. In 
the simulations, H\ cannot be smaller than the typical 
distance between the discrete values of the dipole field, 
which has a lower limit since the sample cannot be very 
large (see Fig. |l|c); however, the sample dimension in our 
simulation allows us to employ a value of H\ of the same 
order of magnitude of that deducible from experiments. 
In evaluating the local field we add to the dipolar con- 
tribution given by Eqs. (@) and (^|) a constant, external, 
applied field which can be varied. In most of the figures 
reporting results of simulations we use the reduced field 
h = H/H, where H = gfi B S/a 3 ~ 80 Oersted, for S = 10 
and a = 13 A. 

The Monte Carlo algorithm employed to simulate the 
dynamics of our systems is an implementation of those 
proposed and discussed elsewhere 

We start from a completely magnetized sample at time 
t = 0, and we evaluate the local field Hi acting on any site 
i of the lattice; we then cycle through the following steps: 

i) We single out those sites where \Hj,\ < Hi; let us de- 
note with no the number of such spins which, according to 
our model, can relax with a probability given by Eq.(|lO|). 
ii) We increment time replacing t by t + At, where At is 
chosen stochastically with probability r] n e- T ' onoAt , i.e. 
we set At = — ln£/ (j]ono), where £ is a generated number 
uniformly distributed in (0,1). in) We randomly choose 
one of those no spins singled out in step i) and flip it 
with probability given by Eq.10. iii a ) If the spin has been 



flipped we update the total magnetization and the fields 
on all sites of the lattice, iv) We come back to step i). 

For any set of the simulation parameters different in- 
dependent runs were made and averaged. Some sample 
runs where also done using the more elementary algorithm 
which uses constant time steps; the results obtained are 
the same, but it takes much longer time, due to the small- 
ness of the time step to be used to get stable results. 



6 Monte-Carlo simulations: results 
6.1 Short time behavior 

6.1.1 Spherical samples 

For a spherical sample, the Prokofiev-Stamp formula ( |l6| ) 
is very well satisfied if the uniform, initial local field H (0) 
is smaller than the width H\, (Fig. ||). For those fields, the 
agreement with ([l6]) is good (except at very short times) 
as long as the magnetization is larger than 80% of its ini- 
tial value. At extremely short times, the magnetization 
decays linearly with t as well known Qj. This can be un- 
derstood because the spins relax independently since the 
average distance between reversed spins is very large and 
the dipole field created by the reversed spins is smaller 
than the width H\. 

Away from the resonance, i.e. for H(0) = H > Hi, 
formula ( [l6| ) is no longer satisfied. However, if the ini- 
tial field is not too far from the resonance, the magne- 
tization becomes a linear function of \/t after a certain 
time ti(H). This phenomenon can presumably be inter- 
preted as follows. For short times, P(r, t; 0) = except 
near the surface, so that dm/dt is small according to Eq. 
(pp[). However, a partly relaxed zone, with P{r,t;Q) ^ 0, 
will progressively invade the whole sample. The invasion is 
total at time T\(H). This time can be expected to become 
very long for big samples. 



6.1.2 Parallelepipedic samples 

We investigated an elongated, parallelepipedic sample whose 
shape (12 x 17 x 36 spins) roughly corresponds to that 
which was experimentally measured ||. The Prokofiev- 
Stamp prediction ( |l9| ) is satisfied for short times as shown 
by figure [|. However, these 'short' times are differently 
short for different external fields. For certain field values, 
the slope \dmjd\ft\ suddenly decreases after a rather short 
times and the magnetization curve crosses other magne- 
tization curves which satisfy the Prokofiev- St amp predic- 
tion on a longer interval. This crossing has not been ex- 
perimentally seen. It is of interest to relate the different 
shapes of the demagnetization curve to the initial internal 
field density p(0; H) shown by Fig. ||a. The curve p(0; H) 
has a sharp maximum, decreases abruptly to on one 
side and much more smoothly on the other side. The field 
values which satisfy the Prokofiev-Stamp prediction for 
a long time correspond to the smooth side. Those which 
correspond to the steep side and to the maximum satisfy 
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the Prokofiev- St amp prediction during a shorter time. We 
have no precise explanation for these observations, but it 
seems to be related to the sharp edge of the distribution 
and to the fact that the spin interested by relaxation are in 
this case mostly concentrated on the surface of the sample. 

6.2 Long time behavior 

6.2.1 Spherical sample 

Figure |^ shows the total magnetization at long times in 
the case of a sphere. It turns out to be well fitted by a 
stretched exponential with j3\ = 0.37; however the good 
agreement with (jfy is not very significant since, as will 
be seen later, the parameter j3i seems to depend on the 
sample shape, and the experiments are mostly done on a 
parallelepipedic sample. No good fit by (|l5|) is possible. 

The distribution p(t; H) of the internal fields, initially 
very sharply peaked, smoothens as shown by Fig. f?|. 

6.2.2 Cubic and parallelepipedic samples: effect of geometry. 

The magnetization of a cube is shown by Fig. |^, together 
with the magnetization of smaller cubes of various sizes 
having the same center as the whole sample. The external 
field H is chosen such that the internal field vanishes at 
the center of the sample at t = 0, i.e. H = in our model. 
The magnetization of the central region is seen to have 
almost completely vanished while the relaxation is still 
very weak at the periphery of the sample. This suggests 
that the latest stage of the relaxation is dominated by the 
motion of the boundary of the relaxed region. 

The cubic shape of the inner regions chosen in Fig. || 
is the simplest choice, but perhaps not the best. At short 
times, the relaxed region is indeed expected to be the vol- 
ume where the internal field vanishes in the saturated sam- 
ple, and this volume is more complicated than a cube. 

The relaxation of the central region is well described 
by the formula 

t = A In 2 \m{t) /s] + B In 3 [m{t) / s] , (22) 

with B/A w 0.1. 

The decay of the total magnetization is very well fit- 
ted, until 0.15 times the saturation value, by a stretched 
exponential with f3\ = 1/4, which is quite far from the 
experimental value @ for elongated samples. A simula- 
tion done in an elongated parallelepipedic sample (Fig. ||) 
yields (3\ = 0.33, which is closer to the experimental value. 

The distribution p(t;H) of the internal fields, initially 
very sharply peaked, smoothens as shown by Fig. [l^ in 
the case of a parallelepiped. 



7 Conclusions 

Our contribution complements the work of Prokofiev & 
Stamp H(| and suggests new experiments. In particular, 
it is shown that the overall magnetization is remarkably 
well fitted by a stretched exponential, with an exponent 
which seems to depend on the sample shape and takes 
values from 1/4 to 0.4. Moreover, it is pointed out that 
the slow relaxation observed in parallelepipedic samples 
is mainly a result of the slow extension of the relaxed 
region, and the central region relaxes much more rapidly. 
The experimental check of this property would be a crucial 
test of the model. It might be feasible using multisquids 
as in the recent experiment of Wernsdorfer et al. J5J. 

Fruitful discussions with D. Gatteschi and R. Sessoli are grate- 
fully acknowledged. 
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6.2.3 Effect of H x 

The effect of Hi is shown by Fig. [II]. The time scale is 
modified as expected, but the dependence of the relax- 
ation time with respect to the external field is not strongly 
modified, and well compares with that observed in the ex- 
periments. 
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Fig. 1. Distribution of local fields at the sites of a simple 
cubic lattice occupied by spins parallel to the z axis, a) Sphere 
of 9134 spins. The wings are due to the surface, as shown in 
picture b). In such picture the clear regions are those where the 
local field is larger than its average value, the dark regions are 
those where it is lower, c) Cube of different sizes: 10 x 10 x 10, 
20 x 20 x 20, and 40 x 40 x 40 spins, starting from below. 
The field is measured in units gfiBS/a s (~ 80 Oersted, setting 
S = 10 and a ~ 13 A, average distance among Fes clusters in 
the actual compound). Note as the size of the sample affects 
the discretization of the field values. 
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Fig. 2. a) Distribution of local fields at the sites of a simple 
cubic lattice occupied by spins parallel to the z axis, for a par- 
allelepiped of 12 x 17 x 36 spins, b) The local fields through this 
sample: below each picture is reported the value of the external 
field to be applied to bring the shaded spins to resonance (the 
local field is given in units g/isS/a 3 ~ 80 Oersted, see text in 
Sect.|). 
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Fig. 3. Distribution of internal fields in a magnetically satu- 
rated sample of cubic shape, a) Spins localized at the Fe-ions 
sites of the actual material, ignoring the actual electronic den- 
sity, b) Single giant spins localized at the sites of a triclinic 
lattice with the unit cell parameters of Fes- c) Single giant 
spins localized at the sites of a cubic lattice with lattice con- 
stant a — 13 A. d) The same distributions given in the previous 
picture, with a) (open squares) and b) (open triangles) shifted 
by 280 and 360 Oersted, respectively. 
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Fig. 4. Magnetization as a function of the square root of the 
time for a spherical sample of radius 13 lattice constants for 
»7(0) = 10 4 and various applied fields h. From the uppermost 
curve: h=0.3, h=0.2, h=0.15, h=0.1, h=0.05, h=0. The re- 
ported data are the average over 10 independent runs. In the 
inset the field distribution is reported. 



Fig. 5. Magnetization as a function of the square root of the 
time for a parallelepipedic sample of 12 x 17 x 36 spins for 
77(0) = 10 4 ; the longest side is along the easy magnetization 
axis. The various curves correspond to the following applied 
fields: a) h=3; b) h=2.6; c) h=3.2; d) h=2.4; e) h=2.2; f) h=3.4; 
g) h=2; h) h=1.6; 1) h=3.6. The reported data are the average 
over 10 independent runs. 
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Fig. 6. Long time relaxation of the magnetization in a spher- 
ical sample of radius r — 13 for 77(0) = 10 4 . Continuous lines: 
magnetization as function of time for external applied field 
h = 0.1 ( uppermost curve, fitted stretched exponential pa- 
rameters: /3i = 0.36, r = 0.036) and h = (fitted stretched 
exponential parameters: (3\ = 0.33, r = 0.037). The dashed 
lines show the fraction of spins which have never flipped, for 
the same values of the applied field and in the same order. 
The thin, steady lines through the magnetization curves are 
the stretched exponential fits. 
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Fig. 7. Evolution of the field distribution in a sphere with no 
applied field; in each picture the value of the corresponding 
magnetization is given. 
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Fig. 8. Magnetization as a function of time for a cubic sample 
of 32 3 spins and in regions of cubic shape having the same 
center as the cube. The external field is chosen such that the 
internal field vanishes at the center of the sample at t = 0, 
Hi = 0.1, r/(0) = 1 (in order to compare with the data for other 
sample shapes reported in the other figures, please note that a 
change of 77(0) entails only a rescaling of time), and the average 
is done over 20 realizations. Upper curve: magnetization in the 
full sample of 32 3 spins. Lower curves: magnetization in sub- 
cubes, concentric to the full one, of the dimension given in the 
legend. 



Fig. 9. Long time relaxation of the magnetization in a par- 
allelepipeds sample 12 x 17 x 36 for 7/(0) = 10 4 in external 
applied field h = 2.4. Continuous line: magnetization as func- 
tion of time; dashed line: fraction of spins which have never 
flipped; the thin, steady line through the magnetization curves 
is the stretched exponential fit with parameters /3i = 0.33, 
t = 0.026. 
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Fig. 10. Evolution of the field distribution in a parallelepiped 
12 x 17 x 36 with external applied field h — 2.4, corresponding 
to a local field hi = —0.704 at the center of the sample; in each 
picture the value of the corresponding magnetization is given. 
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Fig. 11. Short time inverse relaxation time r _1 as a function 
of applied field for a sphere of radius r = 13 (lower pictures) 
and a parallelepiped 12 x 17 x 36 for two different widths of the 
resonance, Hi = 0.1 (left) and Hi = 0.2 (right) and 77(0) = 10 4 . 



